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Abstract. 

We report the creation of ultracold bosonic ^^Na®^Rb Feshbach molecules via 
magneto-association. By ramping the magnetic field across an interspecies Feshbach 
resonance, at least 4000 molecules can be produced out of the near degenerate ultracold 
mixture. Fast loss due to inelastic atom-molecule collisions is observed, which limits 
the pure molecule number, after residual atoms removal, to 1700. The pure molecule 
sample can live for 21.8(8) ms in the optical trap, long enough for future molecular 
spectroscopy studies toward coherently transferring to the singlet ro-vibrational ground 
state, where these molecules are stable against chemical reaction and have a permanent 
electric dipole moment of 3.3 Debye. We have also measured the Feshbach molecule’s 
binding energy near the Feshbach resonance by the oscillating magnetic field method 
and found these molecules have a large closed-channel fraction. 
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The creation and manipulation of ultracold heteronuclear molecules have received 
intensive attentions in recent years due to the versatile and promising potential 
applications [1, 2] of these molecules. With controllable, anisotropic and long range 
dipole-dipole interactions, they could be used in quantum computation [3, 4], quantum 
simulation [5], precision measurement [6, 7] and controlled cold chemistry [8]. So far, 
the most successful scheme for producing ultracold ground-state dipolar molecule is 
by associating ultracold atoms near Feshbach resonances [9, 10] to form weakly-bound 
molecules first, followed by a stimulated Raman adiabatic rapid passage (STIRAP) [11] 
to transfer them to a deeply bound state [12, 13]. This has been successfully applied to 
the ^°K®^Rb system [14], where near degenerate ground-state dipolar fermionic molecules 
are created. However, the chemical reaction 2 KRb—)■ K 2 + Rb 2 is an exoergic process 
which results in a large inelastic loss, severely limiting the trap lifetime of the KRb 
molecular gas [15, 16, 17]. 

Currently, there is a great effort in generalizing the KRb production scheme to other 
heteronuclear alkali dimers. Creation of Feshbach molecules of RbCs [18], LiNa [19] and 
NaK [20] were already reported in 2012, and very recently ground-state RbCs molecules 
were successfully produced [21, 22]. In this work, we focus on the bosonic NaRb 
molecule, which in the absolute ground state is stable against chemical reactions [23] and 
has a permanent electric dipole moment as large as 3.3 Debye [24]. The ground-state 
NaRb molecule can be readily polarized with a moderate electric held. For instance, at 
5 kV/cm the induced dipole moment is already more than 2 Debye. Therefore, it is an 
appealing system for studying the bosonic quantum gas with strong dipolar interactions. 
Recently, the double species Bose-Einstein condensates (BECs) of ^^Na and ®^Rb atoms 
have been produced [25] and their interspecies Feshbach resonances(FRs) [26] were also 
investigated in our group. One of the s-wave resonances between atoms in their lowest 
hyperhne Zeeman state locates conveniently at a magnetic held Bq = 347.7 G with 
a width of A = 4.9 G. Here we report the creation and characterization of ^^Na®^Rb 
Feshbach molecules by magneto-association with this FR. 

We prepare the ultracold mixture of ^^Na and ^^Rb atoms in the same setup for 
producing the double species BEG [25]. Here we only describe it briehy. The initial 
evaporative cooling is performed in a hybrid trap [27] with ^^Rb as the coolant and 
the minority ^^Na is sympathetically cooled. The sympathetic cooling is quite efficient 
thanks to the ideal heteronuclear scattering length of 66.8 Oq [26], with oq the Bohr 
radius. The atoms are then transferred to a crossed optical dipole trap (ODT) where 
the hnal evaporative cooling is performed. In our previous work, the ODT was formed 
by a 1070 nm laser which generates very different trapping potential depths for Na and 
Rb, with a ratio ~ 1/3. Thus when lowering the ODT power, Na always 

evaporates faster than Rb, limiting the double BEG numbers to about 4x10^ for each 
species. Here, we add in another single beam trap created by a 660 nm laser (Gobolt 
Flamenco), which generates a repulsive potential for Rb and an attractive one for Na. 
In this dichromatic ODT, we can tune the relative trap depths almost at will to allow 
the majority Rb atoms to always act as the coolant in the ODT and minimize the loss 
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of Na atoms. With this upgrade, the double BEC atom numbers are increased to more 
than 10^ for both species. 

For the best atom-to-molecule conversion, it is crucial for the two atomic clouds 
to have the maximum overlap in phase space [28], which ideally should be achieved 
with both species condensed. However, we have found empirically in our experiment 
that the largest pure Feshbach molecule sample can only be obtained with a near¬ 
degenerate mixture. This is believed to be a result of the very fast inelastic losses 
caused by collisions with the high density residual atoms inherent from the condensates. 
We also observe that better results can be obtained by turning off the 660 nm beam 
adiabatically at the end of the evaporative cooling. This is understood as a result of the 
smaller differential gravitational sag and thus better spatial overlap between the two 
clouds, as trap oscillation frequencies for the two species are much more similar in the 
pure 1070 nm ODT than in the dichromatic trap. 

During evaporation, both atoms are transferred to their \F = l,mF = 1) hyperhne 
ground state with a radio-frequency adiabatic rapid passage in a 2 G magnetic field. 
The magnetic field is then ramped up to 364 G, far above the 347.7 G FR. Evaporation 
continues at this field and stops near the critical temperature of the Na BEG. We then 
jump the magnetic field down to 350 G, the starting point of the magneto-association. 
After 5 ms holding for the magnetic held to stabilize, we typically have an ultracold 
mixture with 8 x 10^ Na and 1 x 10^ Rb atoms at a temperature of 350 nK. The average 
hnal trap oscillation frequency is measured to be cJ = 27r x 170(148) Hz for Na(Rb). 
The calculated peak number density and phase space density are 4 x 10^^(2.6 x 10^^) 
cm“^ and 1.0(0.84) for Na(Rb), respectively. 




Figure 1. Magneto-association and dissociation, (a) Remaining atoms after the 
magnetic field is swept downward and stopped at various values near the FR. (b) 
Reverse the sweeping direction without holding after the association, an atom recovery 
shows up as a result of molecule dissociation. Due to limited reverse ramp rate, we 
stop the association close to resonance at 347.2 G to shorten the dissociation sweep 
duration. Red solid lines are fit to the hyperbolic tangent function. Error bars are 
from statistics of typically 3 shots and represent one standard deviation. 

To perform magneto-association, we adiabatically sweep the magnetic held toward 
the FR and across it at a hxed rate of 5.2 G/ms. As shown in Fig. 1(a), we stop at 
various held values and then turn the held oh rapidly without holding before measuring 
the number of remaining atoms. Huge losses are observed for both species after the FR 
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Figure 2. Creation of pure Feshbach molecule samples by removing residual atoms 
with a high magnetic field gradient after association, (a) shows the magnetic field 
sequence and (b) shows images of Na and Rb at different stages of the association 
and dissociation procedure with the removal gradient pulse always present (see text for 
details). Color bar denotes optical density. By varying the holding time before or after 
the gradient pulse, we can measure the lifetime of Feshbach molecules with or without 
residual atoms. 


is crossed, amount to about 50% or 4 x 10^ atoms for the minority Na atom. We note 
that these losses are the combined result of Feshbach molecule creation and enhanced 
three-body recombination. To conhrm Feshbach molecules are indeed created, we ramp 
the held back up right after the association with a hxed rate of 3.9 G/ms. During this 
process, any Feshbach molecules will be dissociated to free atoms and show up as an 
atomic number recovery. A typical dissocaition curve is shown in Fig. 1(b). An increase 
in the number of atoms is clearly observed after the held is ramped up across the FR, 
which verifys that Feshbach molecules were created and then dissociated. However, only 
10% or 4000 of the lossed atoms were recovered. We believe that this is a lower bound 
of the Feshbach molecule number as fast atom-molecule collisions during the hnite ramp 
up time can kill these molecules. 

In future molecular spectroscopy studies and quest for absolute ground-state NaRb 
molecules, it is desirable to have a pure sample of Feshbach molecules. To this end, we 
have to remove residual Na and Rb atoms quickly after the association. According to the 
previous coupled-channel modeling [26], Feshbach molecules created with the 347.7 G 
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FR have a nearly zero magnetic dipole moment when the binding energy is more than 
11 MHz which corresponds to magnetic fields ~7.5 G below resonance. As residual 
atoms in |1,1) state are high field seekers, we can thus remove them with a magnetic 
field gradient without affecting molecules. As illustrated in Fig. 2, this is accomplished 
by a 3 ms gradient pulse of 150 G/cm with the Feshbach magnetic field at 335.0 G. 
After this, molecules are dissociated and the magnetic field is turned off before images 
are taken. With the high magnetic gradient applied, we have confirmed that both 
Na and Rb atoms can be removed completely without the association and dissociation 
ramps, or with the association ramp only, as evident by blank images in the first two 
columns of Fig. 2 (b). Signals appeared in the “Association & dissociation” column of 
Fig. 2 (b) which are obtained after completing the full magnetic field sequence, can then 
only come from pure Feshbach molecules created by the association and survived the 
magnetic field gradient. In the optimized condition, we can produce a pure sample 
of 1700 NaRb Feshbach molecules which corresponds to a overall atom-to-molecule 
conversion efficiency (including contributions from both association and inelastic loss) 
of 2%. 



Figure 3. Lifetime of Feshbach molecules. Black solid circles are measured number 
of Feshbach molecules vs. holding time with residual Rb and Na atoms, while red open 
circles are without. Solid curves are exponential fitting for extracting lifetimes. Error 
bars represent one standard deviation from typically 3 shots. All measurements are 
performed at a magnetic field of 335 G, where Feshbach molecules are almost of pure 
closed-channel character (see discussions below). 


We notice that similar overall conversion efficiencies were also reported by 
other groups working on different bosonic heteronuclear molecular species [18, 29]. 
Nevertheless, we want to emphasize that this low molecule number is not due to the 
association step only, which actually could have a much higher conversion efficiency. It 
is also a result of inelastic atom-molecule collisions. Even with the high field gradient, 
we estimate that about 1 ms is still needed for separating molecules from residual atoms 
completely, which is still ample time for inelastic collisions. As shown in Fig. 3, we 
have measured the Feshbach molecule lifetime with and without the residual atoms by 
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varying the holding time before and after the gradient pulse. With residual atoms, the 
lifetime is 1.1(1) ms. In 1 ms, almost 60% of the population will be destroyed. With 
the residual atoms removed, the lifetime increases by a factor of twenty to 21.8(8) ms. 

This lifetime and the current signal to noise ratio are good enough for further 
experiments towards ground state molecules. Faster residual atom removal with light 
blasting [30, 14] should be able to increase the number of pure molecules. Ultimately, 
inelastic losses may be eliminated by hrst loading both atoms into three dimensional 
optical lattices and making a double species Mott insulator before the association 
[30, 31, 32, 33]. Currently, our numbers of atoms in the atomic mixture are also on 
the low side, but we cannot see any fundamental limitations which will prevent us from 
increasing them signihcantly by technical improvements. 




Figure 4. Feshbach molecule binding energy measurement by the oscillating B 
field method. Inset of (a): an oscillating field can transfer atom pairs to a molecular 
state when its frequency matches with the Feshbach molecule’s binding energy which 
is determined by the Feshbach field detuning. When scanning the Feshbach field in the 
presence of a 1 MHz oscillation field, besides the bare resonance at 347.7 G, additional 
loss features show up for both Rb (a) and Na (b) when the resonant condition is met. 
Solid dots are experimental data and solid curves are fitting results (see text). The 
dashed vertical lines indicate the resonance magnetic field found from the fit. 

To further characterize the Feshbach molecule, we measure its binding energy vs. 
the magnetic held detuning with the oscillating held method [34]. The oscillating held 
is produced by a single loop coil driven by a 2 W radio frequency ampliher. It is placed 
coaxially with the Helmholtz coils producing the Feshbach held. As illustrated in the 
inset of Fig. 4(a), when the oscillation frequency matches with the Feshbach molecules’ 
binding energy at the selected magnetic held detuning, transition between the free atom 
and Feshbach molecule states will be induced. This coupling can associate pairs of free 
atoms to Feshbach molecules, which will get lost from the trap quickly due to inelastic 
collisions. 

As the single loop coil has a rather narrow bandwidth, maintaining a constant 
held amplitude over all relevant frequencies is difficult. Similar to reference [18], for 
the binding energy measurement, we instead hx frequency of the oscillating held, and 
scan the Feshbach magnetic held. As shown in Fig. 4(a) for Na and (b) for Rb are 
such a measurement in the presence of a 1 MHz oscillating held. Besides the main 
broad loss feature centering at 347.7 G, another narrower dip shows up on the lower 
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Figure 5. Binding energy vs. magnetic field near (a) 347.7 G and (b) 478.8 G FRs. 
Open circles are experimental data. Red solid lines are fitting results with the square 
well model and dashed lines indicate the thresholds. Error bars on magnetic field are 
from fitting and magnetic field calibration. The inset of (a) shows the closed-channel 
fraction of the Feshbach molecule vs. magnetic field near the 347.7 G FR calculated 
from the fitting parameters. 


field side when the Feshbach molecule’s binding energy matches with the oscillation 
frequency. To extract the resonant magnetic held at this oscillating frequency, the 
apparently asymmetric lineshapes are htted with a Gaussian function convolved with the 
Boltzmann distribution due to the hnite temperature of our sample [35, 36]. Typically, 
the center of the measured lineshape is shifted away from the resonance, as shown by the 
vertical lines in Fig. 4(a) and (b). The same measurements are repeated for oscillation 
frequencies from 100 kHz to 3.5 MHz, beyond which the coupling becomes too weak for 
reliable measurement due to reduced free-bound Franck-Condon factors [37]. 

The results are summarized in Fig. 5(a). We have also performed the same 
characterization for another s-wave FR at 478.8 G [26], as shown in Fig. 5(b). We notice 
immediately that binding energies change quickly with magnetic held detuning for both 
resonances. Within 2 G, the Feshbach molecules are already bound by more than 3 MHz, 
which is a strong evidence of large closed-channel fraction or that these molecules are 
“non-universal” [9] . This hnding is not surprising as the background scattering length is 
a positive and moderate value of 66.8 oq, which indicates that closed channels for both 
FRs are real molecular states. We note that this behavior agrees qualitatively with the 
prediction from the coupled channel calculation [26]. 

To extract quantitative information from the binding energy measurement, we £t 
our data with the square well model developed by Lange et al [38]. Near a Feshbach 
resonance with intermediate width, the binding energy Ef, = can be determined 
from the magnetic dipole moment difference between the open and closed channels 
and the Feshbach coupling strength F as 

1 F/2 

^ a {Ek + 6fi {B - Bf)) ■ 

Here is the wave number. Be is the magnetic held where the bare molecular state 
is tunned to the open-channel threshold and d is the mean scattering length which can 
be calculated from the van der Waals Cq coefficient. The solid curves in Fig. 5 are the 
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Table 1. Parameters obtained from the square well model fitting of the binding 
energy vs. magnetic field data. Results from the coupled channel calculation in [26] 
are also listed as a comparison, /tb is the Bohr magneton. 


coupled channel 

square well 

Bo(G) 

A(G) 

Bo(G) 

A(G) 


347.75 

4.89 

347.64(3) 

5.20(0.27) 

2.66(29) 

478.79 

3.80 

478.83(3) 

4.81(0.27) 

2.52(26) 


fitting results with 6fi, F, and Be as the htting parameters. Here Ce = 1.2946 x 10^ 
cm“^A® is used for calculating a and Obg is taken to be 66.77 Oq [26]. 

The resonance width A = and resonance magnetic held Bq = 

B -then be obtained from combination of these parameters, as summarized 

in Table. 1. While Bq for both resonances agree well with results from the global 
coupled channel modeling, the resonance widths show some discrepancy. We also want 
to mention that binding energies predicted from coupled channel calculations also have 
some disagreement with our measurement. We hope to resolve these problems by adding 
more data points to the Feshbach spectroscopy with improved resolutions later. 

From the htting parameters, the widely used dimensionless resonance strength 
parameter Sres [37] for the two FRs are calculated to be 0.72 and 0.63, respectively, 
which conhrms that their coupling strengths are both in the intermediate regime. We 
have also deduced the closed channel fraction of the Feshbach state. As shown 

in the inset of Fig. 5(a), for the 347.64 G FR, the closed-channel character increases 
rapidly with the magnetic held detuning. At about 346 G, the closed-channel fraction 
is already 50% and reaches over 80% at 340 G. 

In conclusion, we have successfully produced ultracold ^^Na®^Rb Feshbach molecules 
by Feshbach magneto-association. After removing residual atoms, the small pure 
molecular sample lives long enough for further work. These Feshbach molecules have a 
large closed-channel fraction near the Feshbach resonance, advantageous for hnding high 
efficiency STIRAP routes for population transfer. We note that a promising path for 
population transfer via the 2^S’''/l^n excited state admixture was already investigated 
in detail with conventional molecular spectroscopy[39]. We are now in a good starting 
point toward chemically stable ground-state bosonic molecules with a large electric 
dipole moment of 3.3 Debye for investigating ultracold gas with dominating dipolar 
interactions. 
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